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Abstract 

We study a statistical model of random plane partitions. The statistical model has 
interpretations as five-dimensional M = 1 supersymmetric SU{N) Yang-Mills on M'^ x 5^ 
and as Kahler gravity on local SU {N) geometry. At the thermodynamic limit a typical 
plane partition called the limit shape dominates in the statistical model. The limit shape 
is linked with a hyperelliptic curve, which is a five-dimensional version of the SU{N) 
Seiberg-Witten curve. Amoebas and the Ronkin functions play intermediary roles be- 
tween the limit shape and the hyperelliptic curve. In particular, the Ronkin function 
realizes an integration of thermodynamical density of the main diagonal partitions, along 
one-dimensional slice of it and thereby is interpreted as the counting function of gauge in- 
stantons. The radius of can be identified with the inverse temperature of the statistical 
model. The large radius limit of the five-dimensional Yang-Mills is the low temperature 
limit of the statistical model, where the statistical model is frozen to a ground state that 
is associated with the local SU{N) geometry. We also show that the low temperature 
limit corresponds to a certain degeneration of amoebas and the Ronkin functions known 
as tropical geometry. 



*E-mail: maedaShet . phys . sci . osaka-u .ac.jp 
^E-mail: nakatsuOhet .phys . sci . osEika-u. ac . jp 



1 Introduction and summary 



It is an old idea that spacetime may have more than four dimensions, with extra coordinates 
being unobservable at available energies. A first possibility arises in the Kaluza-Klein theory. 
In the Kaluza-Klein approach, gravitation and electromagnetism could be unified in a theory of 
five-dimensional geometry. The extra dimension makes that radius is microscopically small. 
The Kaluza-Klein like approach has always been one of the most intriguing ideas in physics. 

The idea of the Kaluza-Klein theory is utilized in superstring theory |IJ. Superstring theory 
is a candidate for a theory of everything. When superstring theory is completed, all four- 
dimensional theories, such as standard model and general relativity, can be derived from super- 
string theory. In superstring theory, strings which are one-dimensional objects play a central 
role instead of point particles. Each elementary particle corresponds to each vibration mode of 
string. While it moves in a spacetime, a string sweeps out a two-dimensional surface. So the 
motion of a string is given by a map from a two-dimensional surface S, called the worldsheet, 
to a spacetime. In this sense the spacetime is often called the target space. Superstring theory 
makes sense only in ten spacetime dimensions (at least in perturbative treatments). Is it a pity 
that the spacetime dimension is not four? We don't think so. We take the ten-dimensional 
space to be of the form M'^ x X, where is our four-dimensional Minkowski space and X is 
a compact six-dimensional space. From the physical viewpoint, the most interesting candidate 
for X is Calabi-Yau threefolds. We cannot look at X directly, however geometrical natures of 
X turn up as physical properties in M^, such as the gauge symmetry and the matter content. 
Studying superstring theory on M'^ x X, one can find that the internal properties of X lead to 
physical consequences for the observers living in M*^. 

Supersymmetric theories [2! give laboratories to test these ideas. Supersymmetric theories 
are more tractable than ordinary non-supersymmetric theories, and many of their observables 
can be computed exactly. Nevertheless, it turns out that these theories exhibit explicit exam- 
ples of various phenomena in quantum field theories. Then it has become a priority in particle 
physics to understand better the perturbative and non-perturbative dynamics of supersym- 
metric theories. M = 2 supersymmetric gauge theories have been particularly studied at both 
perturbative and non-perturbative levels. The extended supersymmetry dramatically simplifies 
the dynamics of gauge theories. Af = 2 gauge theories also have an interesting interpretation 
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from the perspective of superstring theories. It is called the geometric engineering Accord- 
ing to the geometric engineering, M = 2 supersymmetric gauge theory is realized by the type 
IIA superstring on a certain Calabi-Yau threefold. Some physical properties in the gauge theory 
is derived from geometrical natures of the Calabi-Yau threefold. The gauge symmetry is closely 
related with the ADE singularity in the Calabi-Yau threefold. For example, the SU{N) gauge 
symmetry is realized by an A^v-i-type singularity. The geometry which leads to the SU{N) 
gauge theory is often called the local SU{N) geometry. The geometric engineering connects 
four-dimensional gauge theories with geometries of internal spaces. This gives an example of 
the gauge/gravity correspondence. It goes without saying that the most important principle 
in physics is gauge theory and general relativity. The duality between these two theories is an 
interesting research area in particle physics. AdS / CFT correspondence [H E] is one of the most 
remarkable examples of this duality. 

In the last few years, there has been great progress in both superstring theories and super- 
symmetric gauge theories. The exact low-energy dynamics of A/" = 2 supersymmetric gauge 
theories has been revealed by Seiberg and Witten They have evaluated exact low-energy 
effective actions of A/" = 2 theories using the holomorphy and a version of the electromagnetic 
duality. The low-energy effective theory obtained by them led to numerous achievements in un- 
derstanding of the non-perturbative dynamics of gauge theories. Their derivation is elegant but 
follows in a somewhat indirect way. Recently, Nekrasov and Okounkov have calculated exact 
low-energy effective actions in a more direct way [71 |H] . They have evaluated the path integrals 
of A/" = 2 supersymmetric gauge theories and give exact formulae for both perturbative and 
non-perturbative dynamics. The Seiberg-Witten solutions emerges through statistical models 
of random partitions. Also, Nakajima and Yoshioka derived independently the Seiberg-Witten 
solutions by taking the algebraic geometry viewpoint jH]. 

Recent progress in superstring theories is attributed to understanding of topological strings. 
Topological string is a simplified model of superstring theories, first proposed by Witten [TH] . 
which capture topological information of the target space. In the past ten years or so, it has 
turned out that topological strings have an enormous amount of applications. Their structure 
is complicated enough to relate them to physically interesting theories, yet simple enough to 
be able to obtain exact results. There are two different variations in topological strings, the 
A-model and the B-model. We are interested in the A-model. Topological A-model string 
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amplitudes for a certain class of geometries can be computed using a diagrammatical method, 
called the topological vertex ^] . One of the most interesting applications of topological 
strings is the geometric engineering. The topological vertex allows us to evaluate the topological 
string amplitude for the local SU{N) geometry. The result of topological strings actually 
reproduces Nekrasov's formula for A/" = 2 SU{N) gauge theory [T31II1]. 

The method of the topological vertex yields an unanticipated but very exciting connection 
between the topological string and a statistical model of random plane partitions ^Hl- This 
connection has a surprising interpretation as quantum gravity. Quantum theory of gravity 
is the holy grail of physics. In quantum gravity spacetime undergoes quantum fluctuations, 
which cause wild fluctuations of the geometry and the topology of spacetime. These quantum 
fluctuations make spacetime foamy at short length scales [12] . This idea is very exciting, but we 
haven't made a precise understanding of quantum gravity yet. It is believed that superstring 
theory gives rise to quantum gravity on the target space. The topological string can also 
generate a quantum gravitational theory. It is expected that the classical part of this target 
space field theory is the Kahler gravity The classical theory seems to receive quantum 

deformations caused by string propagation. It is conjectured in [TH] that the statistical model 
of random plane partitions is nothing but the quantum Kahler gravity. The gravitational path 
integral involving fluctuations of geometry and topology on the target space is interpreted as 
the statistical sum in the plane partition model. Namely, each plane partition corresponds to 
a geometry of the target space. Then the plane partition model can give a precise description 
of the target space quantum gravity. 

In this article, we investigate a certain model of random plane partitions and its physical 
applications. Our motivation is to understand the relation between superstring theories and 
gauge theories, and clarify the gauge/gravity correspondence between the gauge theories and 
the internal space gravities. We think that the well-defined statistical model gives a useful tool 
to study these relations. Besides, our statistical model would be a good laboratory for studying 
quantum gravity. 
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(a) 



(b) 



Figure 1: The three-dimensional Young diagram (a) and the corresponding sequence of partitions 
(the two-dimensional Young diagrams) (b). 



A plane partition tt is an array of non-negative integers 

TTll 7ri2 7ri3 ■ ■ ■ 
7r2i 7122 7723 " " " 



(1.1) 

TTai 7r32 7r33 ■ ■ ■ 



satisfying vTjj > tTj+ij and TTjj > nij+i for all i,j > 1. It is identified with the three-dimensional 
Young diagram as depicted in Figure l-(a). The three-dimensional diagram vr is a set of unit 
cubes such that vTjj cubes are stacked vertically on each (i, j)-element of tt. The diagram is 
also regarded as a sequence of partitions '/r(m), where m G Z. A partition is identified with the 
(two-dimensional) Young diagram. See Figure l-(b). 

Among the series of partitions, a partition at the main diagonal m = 0, denoted by vr(0), 
will be called the main diagonal partition of vr and play an central role in our argument. We 
consider the following model of random plane partitions. 

TT 

where q and Q are indeterminates. |7r| and |vr(0)| denote respectively the total numbers of 
cubes and boxes of the corresponding diagrams. The model with Q = 1 is well-known [THj . 



The authors of ^3] investigated the Q = 1 model and proposed a connection between the 
model and topological strings. Our plane partition model contains a new parameter Q. By an 
identification of q and Q with the relevant string theory parameters, the partition function (jl.2|) 
can be converted into topological A-model string amplitude on a certain Calabi-Yau threefold. 
We can also retrieve Nekrasov's formulae for five- dimensional Af = 1 supersymmetric Yang- 
Mills theories from the partition function |2Uj . 

The above model has an interpretation as random partitions or the g-deformation. It can 
be seen by rewriting the partition function as 

A 7r(0)=A 

where the summation over plane partitions in ()1.2j) are divided into two branches. The partitions 
A are thought as the ensemble of the model by summing first over the plane partitions whose 
main diagonal partitions are A. 

It is known that a partition A has an alternative realization in terms of charged 
partitions Pr)}^=i- The charges Pr are subject to the condition J2r=iPr = 0- Among 

partitions, these coming from the charged empty partitions {(0,pr)}^i turn to play special 
roles. They are called A^ cores. The above realization allows us to read the summation over 
partitions in p.3|) by means of A^ charged partitions. Regarding the model as the g-deformed 
random partitions, we will factor the partition function into 

Z{q,Q) = 5^Zg-J^)({M;g,g) ^ Z-^[^)({AM},W 

{■Pr} {AM} 

= 5^ ^5t/(7v)({p.};g,Q). (1.4) 

{Pr} 

In the first line we write the Boltzmann weight for the A^ core, that is read from ()1.3p . by 
^sulN)({Pr}'^1^Q)- second line we have included the summation over partitions A^*"^ 

implicitly in Zsu(N){{Pr}', q, Q)- This factorization turns out useful to find out the gauge theo- 
retical interpretations. The relevant field theory parameters are a^, A and i?, where are the 
vacuum expectation values of the adjoint scalar in the vector multiplet and A is the scale pa- 
rameter of the underlying four- dimensional theory. R is the radius of S*^ in the fifth dimension. 
We identify these parameters with g, Q and pr in ()1.4j) as follows. 

g = e-#^ Q = {RK)\ pr = ar/h. (1.5) 
5 



The parameter h is often identified with the string couphng constant gst- The above identifica- 
tion leads to 

ZsU{N)i{Pr}; g, Q) = ^5dSYM({ar}; A, R, K), (1.6) 

where the RHS is the exact partition function [8^ for five-dimensional A/" = 1 supersymmetric 
SU{N) Yang-Mills with the Chern-Simons term. The five-dimensional theory is living on 

X S"^. The Chern-Simons coupling constant Ccs is quantized to A^. Actually, the identification 
()1.5|) shows that the perturbative part and the instanton part of the exact partition function 
for the gauge theory are given by -2'^^*^^ and X]{a{'-)} ^suIn) respectively. 

The gauge theory prepotential is revealed from the exact partition function by taking the 
semiclassical limit, that is, the ^ — > limit. This corresponds to the thermodynamic limit of 
the statistical model. At the thermodynamic limit, the typical volume of the three-dimensional 
Young diagrams or plane partitions is ~ and the variance of the volume is ~ /i"^. Rescaling 
in all directions by a factor h, the typical three-dimensional Young diagrams approach a smooth 
limit shape which is a two-dimensional surface in the octant. For the Q = 1 model, such a limit 
was considered in [221 '^^^ obtain the limit shape as in Figure |21 




Figure 2: The limit shape of the simple model Z(q, Q = 1) = gl'^L 
In the Seiberg-Witten approach, hyperelhptic curves known as the Seiberg-Witten curves 
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play an important role. Given a gauge group and matter content, the potential of the M = 2 
supersymmetric gauge theory has flat directions, so that the M = 2 theory has a space of 
physically inequivalent vacua. This space is often called the vacuum moduli space. The vacuum 
moduli space can be identified with the moduli space of the Riemann surfaces. For example, 
each vacua of the SU{N) Yang-Mills theory corresponds to a hyperelliptic curve with genus 
A^- 1 [2Sj. 

The connection between the plane partition model and the SU{N) Yang-Mills suggests a 
relation between the limit shape and the hyperelliptic curve. In order to find out the relation, 
we focus our attention on the main diagonal partitions and consider the thermodynamic limit 
of each component that appears in the factorization ()1.4|1 of the plane partition model. We 
introduce a density p{x\\;p) of a charged partition (A,p). To obtain a finite limit shape we 
must scale partitions in a certain manner at the thermodynamic limit. Thereby the density of 
a partition A is scaled to p{u\X). The asymptotic form of the Boltzmann weight of the plane 
partition model is expressed as an energy functional of the scaled density. 

iog(gi^| J2 '?'^') = -^{^[p(-lA)] + 0(^)1 (1.7) 

^ 7r(0)=A ^ 

The minimizer of the energy functional gives the typical shape of the main diagonal partitions 
at the thermodynamic limit. The minimizer that realizes the thermodynamic limit of the 
component must be found from configurations that are expressible in terms of charged 
partitions with the fixed charges. 

The variational problem of the energy functional can be solved by using the standard argu- 
ment The solution pf^^^^ turns to be realized using data of the hyperelliptic curve. 

1 ^ 

where /3r are real positive numbers determined by the charges. This curve is thought to be a 
five-dimensional version of the SU{N) Seiberg-Witten curve and reduces to the Seiberg-Witten 
curve of four- dimensional Af = 2 SU{N) Yang-Mills at the limit R 0. 

The more direct link between the limit shape and the Seiberg-Witten curve can be revealed. 
To this end, we introduce a new object called an amoeba fIE\ 1^ . The amoeba of a Laurent 
polynomial f{x,y) is, by definition, the image in of the zero locus of f{x,y) under the 
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simple mapping Log that takes each coordinate to the logarithm of its modulus. We put 
{u, v) = Log(2;, y) to be coordinates of where amoebas live. The most important tool to 
study amoebas is the remarkable Ronkin function Nf{u,v), which is a convex function over 
M^. It has been verified in P2] that the limit shape of the statistical model with Q = 1 is 
expressed by the Ronkin function of an amoeba. We give pieces of evidence for the relation in 
our statistical model. For the SU{N) gauge theory, we employ the polynomial. 

N 

fsuiN){x, y) = l[{x - (3r) - {RAfiy + y-') . (1.9) 

r=l 

The curve fsuiN){^^^,y) = is the hyperelliptic curve ()1.8|1 . We focus our attention on the 
Ronkin function over the M-axis. We will verify the following relation between the Ronkin 
function and the minimizer. 

/u 
dupf'^''\u) + const. (1.10) 
-oo 

Since the minimizer expresses the gradient of the limit shape at the main diagonal, the 
above relation shows that the Ronkin function of fsu(N) is identical with the limit shape of 
plane partitions at least over the main diagonal. In terms of the internal space geometry, the 
limit shape corresponds to a low-energy solution of the Kahler gravity. On the other hand, 
the Seiberg-Witten curve describes the non-perturbative vacuum structure of the M = 2 gauge 
theory. Therefore the above connection gives the full gauge/gravity correspondence. The low- 
energy geometry contains large quantum fluctuations from the local SU (N) geometry. These 
quantum fluctuations would correspond to the instanton correction in the gauge theory. 

One of the most interesting phenomena in superstring theory is the mirror symmetry P7|l28j. 
The mirror symmetry is a T-duality, in which two very different internal spaces give rise to 
equivalent superstring theory. In terms of topological strings, the mirror symmetry states 
that the topological A-model on a Calabi-Yau threefold is equivalent to the B-model on a 
mirror Calabi-Yau threefold. For the local SU{N) geometry, the mirror Calabi-Yau threefold 
is specified by the hyperelliptic curve fsu{N){x,y) = [SH]- Therefore, the above relation is 
also a manifestation of the mirror symmetry. 

Ronkin's functions appear in the dimer problem of bipartite lattice graphs |^ . A dimer 
is a dumb-bell shaped molecule that occupies two adjacent lattice sites, and the dimer problem 
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is to determine the number of ways of covering the lattice with dimers so that all sites are 
occupied and no two dimers overlap. When dimers are magnetically charged, the free energy 
per area is given by a Ronkin's function at the thermodynamic limit, where the coordinates 
{u,v) represent constant magnetic field B [22] • The connection with quiver gauge theories has 
been argued in ||33j. 

The parameter R is identified with the radius of in the fifth dimension of the gauge 
theories. In the plane partition model, if one keeps RA fixed, this parameter is interpreted 
as the inverse temperature, that is, R = 1/T, where T denotes the temperature. Therefore 
the large radius limit of the gauge theories corresponds to the low temperature limit of the 
statistical models. As the temperature approaches to zero, the statistical models get to freeze 
to the ground states. Each ground state is a specific plane partition that is determined by the 

core at the main diagonal. Crystal is the complement of the ground state in the octant. 
It has an interpretation as the gravitational quantum foam of the local SU{N) geometry 
The crystal also has an interpretation in amoebas. We will show that the low temperature 
limit R oo corresponds to a certain degeneration of the amoeba and the Ronkin function 
known as tropical geometry jSHl EH EZl- For instance, facet of the crystal is described by a 
piecewise linear function which is the degeneration of the Ronkin function ^'^^^^^^(u, v) at the 
limit R oo. 

This article is organized as follows. After introducing amoebas briefly in section |2l we asso- 
ciate amoebas with local Calabi-Yau geometries that geometric engineer supersymmetric gauge 
theories having eight supercharges, and investigate the amoebas and their Ronkin functions. In 
sectioniniwe study the thermodynamic limit of the plane partition model. The relation between 
the limit shapes and the Ronkin functions is shown. In section 0] we describe degenerations 
of the amoebas and the Ronkin functions, and relate the degenerations with crystals that are 
realized at the low temperature limit or the large radius limit. 

2 Amoebas from local Calabi-Yau geometries 

We associate amoebas with local Calabi-Yau geometries that geometric engineer [3] supersym- 
metric gauge theories having eight supercharges. The amoebas live in M^. We investigate the 
Ronkin functions for these amoebas. The Ronkin function is a convex function that is strictly 
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convex over an amoeba and piecewise linear over the amoeba complement. We start this section 
with an introduction to the notion of amoebas |25J and related objects. 



2.1 What is Amoeba ? 

Let f{x,y) = 'Ylii jez^ij^^y^ ^ Laurent polynomial with only a finite number of the flj^'s 
being non-zero. Let V/ be the zero locus of the polynomial in (C*)^, that is, 

Vs = [{x,y)e{C*f\ /(x,2/) = 0}. (2.1) 

This Vf is a punctured Riemann surface. Define the logarithmic map by 

Log : (C*)2 — . IR2 

(2.2) 

(x, y) I — > (log|x|, log |?/| ) . 
The amoeba of / is the image of Vf. 

Af = 'Log{Vf). (2.3) 

The amoeba will typically be a subset of with tentacle-like asymptotes going off to infinity 
and separating the complement '^Af = M.^\Af into some connected components. Every con- 
nected component of the amoeba complement is a convex domain in M^. Simple examples of 
amoebas are shown in Figure El 

One of the main tools to study amoebas is the Ronkin function [26j. We now take {u,v) 
to be Log{x,y), coordinates on where amoebas live. The Ronkin function is defined by the 
integral 

Nf{u,v) = -^[, ——\og\f{x,y)\. (2.4) 
' ' 27ri 2 J\x\=e- X y ^^■'^ ^ ' 

|j/|=e" 

The above integrations are over the preimage of a point {u, v) G under the logarithmic map 
()2.2|) . It is clear that the integrand is singular in the amoeba. But the Ronkin function takes 
real finite values even over there. In fact, Nf is a convex function which is strictly convex 
over Af and linear over each connected component of '^Af. This particularly means that the 
gradient VNf has a definite value over each connected component of ^.4/. As an example, the 
Ronkin function of 1 + x + ?/ is plotted in Figure ^ 
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(a) 



(b) 



Figure 3: The amoebas of (a) 1 + x + y and (b) 5 + x + x ^ — y — y ^■ 

The amoeba and the Ronkin function reflect structure of the Newton polygon of the Laurent 
polynomial. The Newton polygon of / is the minimal convex polygon which contains all points 
G corresponding to monomials in the polynomial. 



The Newton polygons of 1 + x + y and 5 + x + x~^ — y — y ^^ can be found in Figure El 

Over each connected component E of the amoeba complement ^^.4/, the gradient ViVy takes 
value in a integer lattice point of the Newton polygon |26j . 



This gives an upper bound to the number of connected components of '^Af. The connected 
components are always less than the lattice points of the Newton polygon. The symbol 
will be used to denote the connected component over which the gradient VNf equals to («, j). 
For almost all the polynomials considered in this paper, the upper bound is attained, that 
is, the amoeba complement has the same number of the connected components as the lattice 
points of the Newton polygon. For such a polynomial /, we have 



A/ = convex hull of {(i, j) G | a^j 7^ O}. 



(2.5) 



VNf 1^ e Aj nz^ 



(2.6) 




U ^(^.^) • 



(2.7) 



{jj)eA^nz2 
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Figure 4: The Ronkin function of 1 + x + y. 



The Ronkin function is piecewise linear over the amoeba complement "^Af. For each con- 
nected component E of ^Af, let A''^ be the linear extension of Nf\E to M^. When E = -E(jj), 
we have 

Ne^^^^^{u,v) = c^ij) + {{u,v), , (2.8) 

where the symbol ( , ) denotes the standard inner product on M^, that is, {{u, v), = iu+jv. 

If the above (z, j) is a vertex of A/, the constant C(ij) is computed and simply becomes C(ij) = 
^(Inaij). According to |26,, we introduce a piecewise linear function Sf over by 

Sf(u,v) = 'ma.xNE{u,v) . (2.9) 

E 

Taking account of the convexity of the Ronkin function, it follows 

Nf{u,v) > Sf{u,v) (2.10) 
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Figure 5: The Newton polygons of {a) 1 + x + y and {b) 5 + x + x ^ — y — y ^■ 

for any {u,v) G M^. The equality holds over the amoeba complement '^Af. So this Sf is a 
piecewise linear function which approximates the Ronkin function by putting sharp corners in 
it. 

The piecewise linear function Sf allows us to define the spine of the amoeba Af. The spine 
S is the corner set of Sj, that is, the set of {u,v) where Sf{u,v) is not smooth. The spines of 
the amoebas in Figure El are drawn in Figure IHl Note that S G Af. We can regard the spine S 
as the 1-skeleton of the convex which is dual to a certain triangular subdivision of the Newton 
polygon Af. 

2.2 Amoeba and the Ronkin function for U{1) theory 

Geometric engineering [S] dictates that local geometries realize supersymmetric gauge theories. 
The topological vertex countings E] of the topological A-model string partition functions 
on local geometries, which involve from the worldsheet viewpoint sums over holomorphic maps 
to the target spaces, support jlSl EI the idea. The local geometries are noncompact toric 
Calabi-Yau threefolds. If a space is toric, many basic and essential characteristics of the space 
are neatly coded and easily deciphered from analysis of the corresponding lattices. As for the 
local geometries, the data of the space are retrieved from a two-dimensional polygon. We will 
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(a) (b) 

Figure 6: The spines of the amoebas in Figure\^ 

regard the polygon as a Newton polygon. Thereby we will define the polynomial for each local 
geometry. 

We first examine an amoeba for the local geometry relevant to the abelian gauge theory. In 
terms of topological strings, this arises from the topological A-model string on 0(B0{—2) P^. 
The toric polygon of this geometry is drawn in Figure [3 

We employ the following Laurent polynomial. 

fuii){x, y)=x-(3-RA{y + y-'). (2.11) 

The parameters (3, R and A in the above are assumed to be real positive numbers. In particular, 
we consider the case of (3 > 2RA. In the field theory language, R is the radius of in the fifth 
dimension and A denotes the scale parameter of the underlying four- dimensional theory. It is 
clear that the Newton polygon of this fu{i) is the polygon in Figure 

We will rescale amoebas by 1/R. This is achieved by using the following rescaled coordinates 
{u, v) rather than the original ones. 

1 11 
{u,v) = — Log{x,y) = ( — log|x|, — \og\y\). (2.12) 

To accord with this, we also rescale the Ronkin functions by 1/R. In the rest of this paper. 
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(0,1) 




(0,-1) 



Figure 7: The Newton polygon for U(l) theory. 

we describe amoebas by using the above coordinates, and call the rescaled Ronkin functions 
simply as the Ronkin functions. 

Let us now look at the amoeba of fu{i)- The zero locus V/^^^j in (C*)^ is a four-punctured 
sphere. If one regards C* an infinite cylinder, it is a double covering of the cylinder (Figure IH}. 
The amoeba Af^^^^ is the image Log(V/^^^j). Thanks to the condition /3 > 2RA, this becomes 
a subset of described by the following inequalities. 

\P -2RAcoshRv\ < e^"" < p + 2RAcoshRv . (2.13) 

The amoeba is depicted in Figure IHl The amoeba spreads four tentacles. These tentacles 
asymptote to the following straight lines and extend to the infinities. 

They separate the amoeba complement "^Af^^^^ into four connected components. The number of 
connected components agrees with the number of integer lattice points of the Newton polygon. 
Taking account of the above rescaling, the Ronkin function A^/j^^ij is given by 

11/" dx dy 
= ^-(^ j\.\=e^^ —— log I Ml) (a;,?/) I • (2.15) 

Regardless of the simple appearance, it is serious to carry out the above integrations. However, 
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Figure 8: The zero locus V/j^|.^j . 

the gradient VA^'/^^^j becomes tractable since we have 

^^fum 1 f dy f dx 



du (27r2)2 J\y\^^Rv y x- P- RA{y + y ^) ' 

^^Mi) _ 1 /" f d{y + y-^) 



dv {2Txi)'^ Jyj.\^^Ru X J \y\^^Rv RK{y + y ^) — x + jd 

In the double integral ()2.16|) . we can find the following residue integral. 

If dx I 1 if \f3 + RK{y + y-^)\<e^^ 



(2.16) 
(2.17) 



2m X- 13- RA{y + y~^) 



(2.18) 



otherwise. 



The integral takes value 1 or 0. If one puts \y\ = e^^ , the value depends on whether 9y = argy 
satisfies a certain condition or not, as seen from ()2.18|1 . This means that the residue integral is 
a truth function of the condition on 6y. The integration in ()2.16p gives a simple integration 
over By. We thus interpret d^Nf^^^^ as an integration of the truth function over 6y. The similar 
interpretation is also possible for d^Nf^^^^ by using ()2.17j) . 

Over each connected component of the amoeba complement, the gradient VA^/^^^j takes 
value in an integer lattice point of the Newton polygon. It may be helpful to see how this 
happens by using the above interpretation. As an example, let us consider the connected 
component E that is described by the inequality e^" > j3 + 2i?A coshi?t;. We have |a;| > 
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V 




Figure 9: The amoeba of fu{i)- 



1/3 + RA{y + y'^)\ for W {x,y) G Log^^(_E). This means that the residue integral (j2.18p equals 
to 1 and we obtain duNj^^^^ = 1 on this connected component. Computation of d^Nf^^^^ goes 
as well. Thereby we obtain VA^'j^^^j = (1,0). The connected component E is -E(i,o)- The other 
three connected components are -E(o,o) and -E(o,±i) ; where the gradient VA^/j^^^^ takes a on each 
Ea- Over the amoeba complement the Ronkin function becomes the following piecewise linear 
function. 



max 

a=(l,0),(0,±l),(0,0) 



where the constants 

C(1,0) = : 



(2.19) 



C(o,±i) = fog , 



(2.20) 
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We examine the Ronkin function over the amoeba. In particular, we focus on the section of 
the amoeba along the u-axis. In the next section we will see that the Ronkin functions over the 
u-axis relate with instanton countings of gauge theories. The section of the amoeba becomes 
a segment / = [-^ log/5~, log/?"*"], where we put = /3 ± 2i?A. Thanks to the condition 
j3 > 2RA, this segment has a finite length. It is convenient to call the segment the band and 
the complement of it the gap. 



Band 
Gap 



\I- 



(2.21) 
(2.22) 



The residue integral ()2.18|) over the gap takes the definite values irrespective of 6y. It always 
gives for M < log/?" and 1 for -u > log/?"*". On the other hand, when u is in the band, the 
residue integral takes value 1 only if 9y satisfies the inequahty. 



(2.23) 



Otherwise it gives 0. Therefore we obtain 



dN 



fu(i) 



du 



1 

2^ 



de„ ■ 1 



cos Oy<- 



1 + — arccos 

TT 



2RA 



(2.24) 



where the branch of the arccosine is fixed by choosing arccos(O) = — 7r/2. To summarize, the 
gradient of the Ronkin function along the w-axis becomes 

u <j^\og/3- 



dN 



du 







v=0 



1 + ;^ arccos ^- 



2RA 



u e I 
Uog/3+ <u. 



(2.25) 



R 



By using this, we plot the Ronkin function over the u-axis in Figure fTUl For this particular case, 
our computation is not limited on the M-axis and can be done over M^. We plot the Ronkin 
function over in Figure ^2 

The expression (I2.25|) is arranged by using the terminology of complex geometry. The zero 
locus Vfjj^^^ is a complex curve in (C*)^. It can be written as 

^1 _ I ,.Rz 



(2.26) 




u 



Figure 10: The Ronkin function of fu{i) over the u-axis. 



where z = u + i6/R is a cylindrical coordinate of C*. The curve is a double cover of the C*. 
The branch points locate at two ends of the band /. The holomorphic function y has a cut 
along the band on the Riemann sheet and takes values at the unit circle there. In particular, 
it is ±1 at the branch points. Using y, we can rewrite ()2.25|) as follows. 



where we choose arg?/(M — iO) so that it increases along the u-axis from — vr to 0. 

2.3 Amoeba and the Ronkin function for SU{N) theory 

Five-dimensional A/" = 1 supersymmetric SU{N) Yang-Mills theory is realized by a local SU{N) 
geometry. The geometry is an ALE space with Ajv-i singularity nontrivially fibered over P^. 
Fibration of the space reflects the Chern-Simons term of this five-dimensional theory. The 
polygon of this geometry is set out in Figure El 

By understanding the above polygon as the Newton polygon, we employ the following 
Laurent polynomial. 




(2.27) 



fsuiN)ix,y) = QNix) - iRAf{y + y-') , 



(2.28) 
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Figure 11: The Ronkin function of fi[{iy 



where Qn{x) denotes a monic polynomial of degree with real coefficients. Let Pi, . . . , Pn be 
roots of Qn{x). 

Q^(x) = x"" + hx^-' + ■ ■ ■ + bN (2.29) 

N 
r=l 

All the roots are assumed to be real positive numbers. They are arranged in numerical order 
< Pi < ■ ■ ■ < Pn. We also make A very small so that Qn{x) ± 2(i?A)^ have A^ distinct real 
positive roots. 

We observe the amoeba of fsu{N)- The zero locus V/^^f^, in (C*)^ is a four-punctured 
surface. It is a hyperelliptic curve of genus A^ — 1 (Figure HHI). The amoeba ^/^^f^j is the 
image Log(V}-g^^j^j) . The hyperelliptic involution y ^ y^^ induces the involution v —v on 
the amoeba. Shape of the amoeba depends sharply on the parameters A and Pr- When the 
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(0,1) 




(0,-1) 



Figure 12: The Newton polygon for SU{N) theory. 

parameters are set as above, the amoeba spreads its four tentacles and these asymptotes are 

v = ±Nu, v = ±5n, (2.31) 
where Sn is a positive constant given by 

The amoeba complement ^^/g^jjv) 1^^^ A^ + 3 connected components that correspond exactly to 
integer lattice points of the Newton polygon A /sj^^^j . Four of them are unbounded components 
lying between the tentacles. All the other components are bounded components enclosed by 
the amoeba. That is to say, there are — 1 holes in the amoeba. An example of the amoeba 
is drawn in Figure EI 

Consider the section of the amoeba along the w-axis. Let us write 

N 

Qn{x) - 2(i?A)^cos^j, = " Pr{Oy)) . (2.33) 

r=l 

All the l3r{0y)^s become real positive numbers since A is made so small. They are distinct and 
are arranged in numerical order < Pi{Oy) < ■ ■ ■ < (3N{0y)- Bands of the amoeba, which are the 
longitudinal section of the amoeba along the w-axis, consist of A^ segments with finite length. 

^ r 1 1 

Bands : \_\lr, = < - log Pr{Oy) G M : < < 27r ^ . (2.34) 

r=l ^ ' 
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Figure 13: The zero locus Vf^^^^^y 

Similarly to the U{1) case, we call the complement of the bands the gap. Two ends of the r-th 
band are ^\og 13^ and ^\og[3~, where we put [3^ = /?r-(0) and P~ = /^^(Tr). The r-th band 
is a segment between log P~ and log P:^ on the w-axis. Whether P~ < (3^ or /3+ < /?~ 
depends on and r. 



(3~ < /3+ if Nr =even 

p+ < if Nr =odd . 
The Ronkin function Nf^^^^^^ is given by 

11 f dx dy 



(2.35) 



R {27ii) 



\y\=e' 



X y 



^og\ fsu(N)ix,y)\ 



(2.36) 



We will concentrate our attention on the Ronkin function over the w-axis. The gradient along 
the u-axis can be written as 

N 

I / rill I rlT 

(2.37) 



dN 



fsU(N) 



du 



^ 1 



dy 



(271-0^ J\y\=\ y J\x\=eR^ X - Pr{Oy) ' 



There are N poles in the above integrand. Each band holds exactly one pole and every pole 
depends on 6y. It is obvious that every summand in (|2.37|) takes the same form as in the 
previous U{1) case. 

Consider the r-th summand in ()2.37p . The residue integral vanishes unless u satisfies the 
inequality e^" > [3r{0y). When u is off the r-th band J^, this condition is irrespective of 6y. We 
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u 



Figure 14: The amoeba of fsu{2)- 



have 



dy 



dx 



if u < mm{ ^log f3+, ^\og f3- ) , 
J\y\=i y Ji,i=,n.x-Pr{ey) 1 ^ xi>max(llog/3+,^log/3-). 



(2.38) 



On the other hand, when u is on the r-th band J^., we have to read the inequahty e > (3r{0y) 
as a condition on 9y. It should be noted that this reading depends on whether j3~ < or 

< (3~ . We first consider the case of (3~ < The function Qn{x) increases over 
Thus the inequahty can be translated to QN{e^^) > QN{Pr{(^y)) — ^{RA)^ cos 9y over the band 
Ir- So we obtain the following condition for 9y. 



(2.39) 
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Therefore the r-th summand of this case becomes 

dy f dx 1 



de,, ■ 1 



= iV-r + l + -arccos||^LJ|. (2,40) 

where the branch of the arccosine is fixed by choosing arccos(O) = (r — — tt. We turn to 
consider the case of /?+ < In this case, Qn{x) decreases over [Pr^Pr]- ^'^ inequahty 
means Qn{€.^^) < QN{Pr{(^y)) = 2(i?A)^ cos^^y over the band 7^. The condition for 9y becomes 



-».>|g;). (2,41) 



The r-th summand of this case leads to the same expression as ()2.40|) . 

By summing up the above ingredients, we obtain the following expression for the gradient 
of the Ronkin function along the w-axis. 



dNf, 



U(N) 



du 



v=Q 



if max(-^ log/?^) <u < niin(-i log/?^!'^^) 
(r = 0,l,--- ,iV), 

(2.42) 



N + ^ arccos | 



QAr(e«") \ if Ue Is 
2{RA) 



N 



1,2,--- ,iV). 



We put = — oo and P^_^_i = cxd in the above. The graph of duNf^,^^^^\^_^ is depicted in 
Figure IT31 Integrating ()2.42|) . we acquire the Ronkin function over the w-axis. The Ronkin 
function is constant over the connected component -E{o,o) of the amoeba complement and its 
height there is C(o,o) = Sn + log(-RA)^. This gives the boundary condition for the integration. 
See Figure Uni 

We also arrange the expression ()2.42p in terms of complex geometry. The zero locus Vf^^^^^^ 
in (C*)^ is written as 

(2.43) 

where z denotes the cylindrical coordinate of C*. This hyperelliptic curve is a double cover 
of the C* with 2N branch points (Figure E))- The branch points locate at the ends of the 
bands ()2.34|) . The holomorphic function y has cuts along the bands on the Riemann sheet 
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Figure 15: The gradient of the Ronkin function of fsu{3)- 

and takes values at the unit circle over the bands. In particular, it is ±1 at the branch points. 
Using y, we can rewrite ()2.42|) as follows. 



N) 



du 



= ^(n + — \ogy{u-iO)), (2.44) 
.=0 ^ ^ 

where we choose a.Tgy{u — iO) so that it increases along the w-axis from —Nn to 0. 

Note added : It was shown in jSHI (see also |22]) that expression of the gradient like (j2.44p 
holds commonly for the Ronkin functions of Harnack curves, where a curve is Harnack if and 
only if the map from the curve to its amoeba is two-to-one over the amoeba interior. This was 
pointed out by A. Okounkov after the original version of this paper was submitted to e-print 
archives. The authors are grateful to him. 
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3 Ronkin's functions and bd SUSY gauge theories 



By identifying the parameters of the Laurent polynomials with the gauge theory order pa- 
rameters in suitable manners, their Ronkin functions relate with instanton countings of five- 
dimensional supersymmetric gauge theories. The so-called Nekrasov formulae for the exact 
gauge theory prepotentials have a description in terms of random plane partitions. The Ronkin 
functions over the w-axis turn to be identified with integrations of scaled densities of the main 
diagonal partitions in the statistical models at the thermodynamic limit. Thus the Ronkin 
functions are the counting functions. 

3.1 A model of random plane partitions 

A partition A = (Ai, A2, ■ ■ ■ ) is a sequence of non-negative integers satisfying Aj > Aj+i for all 
i > 1. Partitions are often identified with the Young diagrams as seen in Figure IT71 The size is 




Xl(X) X2{X) xi^x){X) 

Figure 17: The Young diagram and the Maya diagram 

defined by | A| = X]i>i which is the total number of boxes of the diagram. A plane partition vr 
is an array of non-negative integers satisfying iTij > Tti+ij and TCij > iXij+i for all i,j > 1. Plane 
partitions are identified with the three-dimensional Young diagrams. The three-dimensional 
diagram tt is a set of unit cubes such that ttij cubes are stacked vertically on each (i, j)-element 
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Figure 18: The 3d Young diagram and the diagonal slices. 

of the diagram. Diagonal slices of vr become partitions. Let 7r(m) denote the partition along the 
m-th diagonal slice. In particular, 7r(0) is the main diagonal partition. This series of partitions 
satisfies the condition 

< 7r(-2) -< 7r(-l) -< 7r(0) y 7r(l) y tt{2) > , (3.1) 

where fi >~ u means the interlace relation between two partitions fi and u. 

IX y V l^l > > H2 > > IJ'3 > ■ ■ ■ ■ (3.2) 

We first consider a statistical model of plane partitions defined by the following partition 
function. 

Ziq,Q) = $^gl-lgW°)l, (3.3) 

where q and Q are indeterminate. The Boltzmann weight consists of two parts. The first con- 
tribution comes from the energy of plane partitions, and the second contribution is a chemical 
potential for the main diagonal partitions. The condition (j3.ip suggests that plane partitions 
are certain evolutions of partitions by the discretized time m. This leads to a hamiltonian for- 
mulation for the model. In particular, the transfer matrix approach developed in [22] express 
the partition function (j3.3|) in terms of two-dimensional conformal field theory. 
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We can interpret the random plane partitions as a g-deformation of random partitions. It 
may be seen by rewriting ()3.3|) as 

^(g,Q) = E^'"( E '?'^')- (3-4) 

A 7r(0)=A 

Partitions A in the above are thought of as the ensemble of the model after summing over plane 
partitions that main diagonal partitions are A. By using the transfer matrix approach [22j, the 
above factorization yields 

Z{q,Q) = Y,Q^'^^x{qKqK---?, (3.5) 

A 

1 3 

where s^{q2, ,■ ■ ■) is the Schur function s^(xi, X2, ■ ■ ■ ) of infinitely many variables specialized 
at Xi = 1^1. 

The statistical model gives rise to a description of five-dimensional A/" = 1 supersymmetric 
Yang-Mills theories To contact with the SU{N) gauge theory, we need to identify the 

indeterminates q and Q with the following field theory variables. 

q = e-#^ Q = {RAf, (3.6) 

where R is the radius of S*^ in the fifth dimension, and A is the lambda parameter of the 
underlying four- dimensional field theory. 

3.2 Thermodynamic limit and variational problem 

The thermodynamic limit is achieved by letting h 0. By summing up partitions A in ()3.5|) . 
we obtain 

Z{q,Q) = l[{l-Qq-r\ (3.7) 

n=l 

The average numbers of cubes and boxes of plane partitions and the main diagonal partitions 
can be computed from ()3.7|) . By using the identification ()3.6|) . the mean values near the 
thermodynamic limit become respectively of orders and h~^. In particular, we have 
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This implies that, as h goes to zero, a plane partition vr that dominates is a plane partition of 
order h~^. Similarly, its main diagonal partition 7r(0) or A becomes a partition of order h~^. To 
realize the thermodynamic limit, it becomes necessary to rescale plane partitions and partitions 
relevantly. 

3.2.1 Scaling partitions 

We provide a description of the scaling of partitions at the thermodynamic limit. Use of the 
Maya diagram becomes helpful. For a partition A, the Maya diagram K{X) C Z is a set of the 
integers Xi{X) = — Aj + i, where i G Z>i. An example of the Maya diagram can be found in 
Figure El For the later convenience, we associate partitions with charges. We denote such a 
charged partition by (A,p), where p G Z is the charge. The Maya diagram K{X,p) is a set of 
the integers Xj(A) +p, that is the parallel transport of K{X) by p along the x-axis. 

For a charged partition (A,p) or the Maya diagram K{X,p), we introduce the density as 



Pbare\'^\ 



|A;p) = 5^5(x-x,(A)-p). (3.9) 

1=1 

The above is not sensitive to the charge since p can be absorbed into the shift of x. We can 
modify ()3.9|) as 

+ 00 +00 

Preg{x\X;p) = ^6{x - Xi{X) - p) - ^6{x - i) , (3.10) 

i=l i=l 

where the subtraction is prescribed so that it satisfies the normal-ordering condition. 

hoo 

dxpreg{,x\X]p) = -p. (3.11) 

X) 

Let A be a partition of order and p be of order h~^. We may think of Xj(A) as quantities 
of order h~^. We regard that elements of such a partition could be suffixed by s G M>o rather 
than i G Z>i, and that these two kinds of indices relate with each other by i = s/h when h is 
nearly zero. As h —>■ 0, Xj(A) are rescaled to a certain function u{s\X), which takes value in R. 

x,iX) = juis\X) + Oih") . (3.12) 



As for the charge we rescale p to a G M by 



p = —a . (3.13) 

/ L 
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It should be noticed that a;i(A), a;2(A), ■ ■ ■ is a strictly increasing series which satisfies the 
conditions that Xi{X) < i for Vi and that Xi{X) = i for i > /(A), where /(A) is the length of 
A. The function u{s\X) becomes a strictly increasing function, and satisfies the conditions that 
u{s\X) < s/N for Vs and that u{s\X) = s/N for s > (, where C is a finite constant obtained 
from /(A). The inverse of u{s\X) exists and is denoted by s{u\X). It becomes a nondecreasing 
function over M, and satisfies the conditions that s{u\X) = Nu for -u > ^ and 5(^1 A) = for 
u < —C,, where ^ is a certain positive constant which one may take ^ = (/N. 

As h ^ 0, we also scale the density ()3.9|1 by using ()3.12|1 and ()3.13|1 . In particular, taking 
account of ()3.12p . we rescale x to u EM.hy x = Nu/h. Then the scaling limit is read as 



lim pbare{x = Nu/k |A; p = Na/h) = p{u — a|A) , (3-14) 



where 



1 ds{u\X) 

pHA) = -^— . (3.15) 

The behavior of the function s{u\X) implies that the scaled density p{u\X) takes values in 
[0, 1]. It also implies that dp{u\X) / du has a compact support in R and satisfies 

= (3.16) 



oo 



du 



Similarly to (|3.14p . the regularized density preg{x\X;p) is scaled to p{u — a\X) — p(m|0). The 
scaled density of the empty partition is p('u|0) = d{u)^ where 9(u) denotes the step function, 
that is, 9{u) = 1 for u > and for u < 0. Applying a partial integration and combining with 
()3.16p . we can translate the condition ()3.11|) as follows. 

^~,™*!m.O. (3.17) 

The scaled density p('u|A) relates with the shape of the Young diagram. A piecewise linear 
curve that traces the upper side of the Young diagram positioned as in Figure El and extends 
over M by |a;|, is called profile of A. The scaling limit of the profile becomes the graph of a 
certain function P{u\X). This function relates with the scaled density by 

^^ = -l + 2pHA). (3.18) 

du 
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3.2.2 The variational problem 

In the description ()3.5|) . each partition A has the Bohzmann weight Q^'^^Sx{q^, q^,- ■ - Y, where 
the parameters are identified with q = exp{—Rh/N) and Q = (RA)"^. When partitions are 
scaled according to ()3.12p . the asymptotic form of their weight is computed in [2^1 by using the 
product formula of the (specialized) Schur function and is expressed as a functional of the 
scaled density. It can be read 

logQi^i.A(gigi---)' = -^{^[p(-|A)]+0(n)}, (3.19) 

where 



E[p{-\\)] = j dudu p{u\\){l - p{u\\))\ogi 



( sinh ^{u — u) 



-oo<«<«<+oo 2 

+iV'^r^.^W). (3.20) 

y_oo 6 du 

Near the thermodynamic limit, partitions of order h~'^ dominate in the statistical model. 
Their Boltzmann weights are measured by using the energy functional E as e~^l^ . Since the 
thermodynamic limit is achieved by letting — > 0, this means that it is realized by a classical 
configuration that minimizes the energy functional. We are thus led to consider a variational 
problem of the energy functional and to find the minimizer as the stationary configuration. 

The minimizer of the energy functional ()3.2()j) should be found from among certain admis- 
sible configurations. The admissible configurations are thought of functions p{u) that could be 
obtained as the scaled densities of partitions. They need to take values in [0, 1]. Furthermore, 
taking account of ()3.16j) and ()3.17|1 . we require that dp{u)/du has a compact support in M and 
satisfies the conditions 

'du^-P^ = 1, (3.21) 
du 

duu^ = 0. (3.22) 
du 

To argue the variational problem, it is convenient to rewrite the energy functional in the 
following form by partial integrations. 

2 Ju^u du du 

. f^^^ Ru'dpiu) 



oo 
oo 



oo 
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where the function 7(m; A; R) has been introduced by the conditions 

9^7, . s /sinh-^Mx 

^(«;A;i?) = log(^^), (3.24) 

7(0; A; i?) = ^(0;A;i?) =0. (3.25) 

It is clear from the expression ()3.23j) that variations p — p + 6p oi the energy functional 
lead to the stationary equation 

The integration symbol in the above means principle part. By differentiating Eq. ()3.26|) twice 
with respect to u, we get 

PP du coth -{u- u)^^ = -1 on supp A . (3.27) 
2 du du 

The variational problem is now stated to find out a solution p^(u) of the integral equation ()3.27|) 

from among the admissible configurations. 

The standard argument [23] allows us to reformulate this as a problem to find out a certain 

analytic function $(^). It may be helpful to consider an analytic function that is realized as 

the integral transform of a continuous function %(«) over the real axis. 



/ + 00 D 
du coth -{z - u) x{u) . (3.28) 



The above function is regarded as an analytic function on C* since it is periodic with the period 
27Ti/R. On the real axis, when one approaches from the upper or the lower half plane, it takes 

<l>^{u ± iO) = PP du coth -{u-u)x{u) T-^xiu)- (3.29) 
J-oo 2 R 

Observing the above formula by putting x = dp/du, we deduce the variational problem as 

follows: Let ^{z) be an analytic function that is periodic with the period 27Ti/R and behave at 

the infinities as 

<^{z) — > ±1 as ^z^±oo. (3.30) 

We further suppose that $(z) has a cut on a bounded region I along the real axis and satisfies 

sR$(m) = -1 on /, 

^ ^ (3.31) 
<^^{u) =0 on M\/, 
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Figure 19: The cut I and the contour C on the Riemann sheet. 

and 

— I dz z^z) = 0, (3.32) 

where C denotes a contour which encircles the above / anticlockwise (Figure Then the 
solution of the variational problem is given by the imaginary part of $ along the real axis 
as follows. 

«cjt(„±,o). (3.33) 

du 2tt 

In particular, we have supp {dpi,/du) = I. 

We note that the conditions (fOT^ and (pO^ are translations of ^JB>, (f^TTfjl and 

()3.17p respectively. 

3.3 U{1) gauge theory and the Ronkin function 

We consider the case of the U{1) gauge theory. We put = 1. The partition function ()3.3|) 
can be identified with the partition function for four- dimensional M = 2 supersymmetric U{1) 
gauge theory on noncommutative R"^. 

To find out a solution of the variational problem posed at the end of the previous subsection, 
we test an analytic function of the following form. 
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where 

l/ + ^/"^ = ;^(e^^-/5)- (3.35) 

Here /? is a real parameter and is assumed j3 > 2RA. 

It is easy to see that the above function satisfies the condition ()3.30|) . Along the real axis, 
it has a single cut on [-^log(/9 — 2RA), log(/5 + 2_RA)]. We can also see that it satisfies 
the condition ()3.3H) . We examine the condition ()3.32j) . Let C be the contour that encircles 
the above cut anticlockwise on the 2;-plane. We will evaluate the contour integral in ()3.32|) as 
follows. 

(p dz z^{z) = —5— (f zdlogy 



^— [ ^\og{RA{y + y-')+p}. (3.36) 

T'-ti J\y\=i y 



2711 

The last integration can be computed by applying the classical Jensen formula in complex 
analysis and gives rise to 



— £ dz zHz) = -- log . (3.37) 



The vanishing of the contour integration fixes f3 as 

(3 = l + {RAf. (3.38) 

We thus obtain the solution p^^^\u) by letting (3 = 1 + (_RA)^ in ()3.35|1 . Let us summarize 
this as follows. 

dp^^^\ , 1 ^ dlogy , 

-^{u) = T-^ —r^i^ ± ^0 , 3.39 

du TT dz 

where 

y + y-' = ^{e^' - 1 - (RAr) . (3.40) 

We note that the above solution is precisely that obtained in |39j by the WKB analysis of the 
fermion wave function. 
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We specialize the Laurent polynomial fu{i), which is given in ()2.1H) . by choosing j3 = 
1 + (i?A)^. The above solution can be expressed by using the Ronkin function of this fu{i)- By 
comparing (j3.39p with (j2.27p . we find 

pUW(u) = ^^'^^^ ^ ■ (3.41) 
ou 

This means that the Ronkin function over the w-axis is an integration of the scaled density of 
the main diagonal partitions that realizes the thermodynamic limit. 

/u 
du p^^^\u) + const. (3.42) 
-oo 

Let P^^^\u) be the limit shape of the Young diagrams obtained from p^^^\u) by using the 
relation ()3.18j) . We can write ()3.42j) as follows. 

0) = i (Pf + + const. (3.43) 

3.4 Random plane partitions and SU{N) gauge theory 

There is a bijective correspondence between charged partitions {p,p) and charged partitions 
{{\^^\ Pr)}^=i- This arises from a division algorithm for the Young diagrams analogous to that 
for integers. The correspondence is neatly described by using the Maya diagrams as follows. 

N 

K{p,p) = \J ir{K{\^'\Vr)) . (3.44) 

r=l 

where '-^ 7j denotes an injective map given by ir{m) = Nm + r — N, and we have 

Lr{K{\^'^\pr)) = N K{X^'^\pr) +r — N. It can be seen that the charges are conserved in (j3.44|) . 

N 

p = J2Pr- (3.45) 

r=l 

In terms of the Young diagrams the correspondence ()3.44|) is seen as Figure EOl 

By using the above correspondence, the factorization ()3.4|) or ()3.5p is also expressed in 
terms of N charged partitions {{X^^\ pr)}^=i- Since the factorization was done by using neutral 
partitions, the charge conservation says that the charges Pr satisfy the SU{N) condition 

N 

J2Pr = 0. (3.46) 



r=l 
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Figure 20: Partition corresponding to three charged partitions {Xi,Pi) = ((2, 1), —7), {X2,P2) = 
((3,2),— 1) and (X^jPs) = ((3, 2), 8) by Three clusters of bold boxes correspond respec- 

tively to the three partitions. 

After including the summation over partitions A^''^ implicitly in each component of the factor- 
ization, let us factor the partition function into 

Z{q,Q) = J2Zsu(N){{Pr};q,Q). (3.47) 

{Pr} 

The above factorization becomes a bridge between the random plane partitions and five- 
dimensional Af = 1 supersymmetric SU{N) Yang-Mills. Let ar be the vacuum expectation 
values of the adjoint scalar in the vector multiplet. We identify the parameters q, Q and Pr 
with the gauge theory parameters as follows. 

g = e-^^ Q = {RAY, pr = ar/h. (3.48) 

It is shown ^21]] that the above identification leads to 

ZsuiN){{Pr}; q, Q) = ^5dSYM({ar}; A, R, K) , (3.49) 

where the RHS is the exact partition function [8^ for five-dimensional A/" = 1 supersymmetric 
SU{N) Yang-Mills plus the Chern-Simons term having the coupling constant equal to A^. The 
five-dimensional theory is living on x S*^, where the radius of in the fifth dimension is 
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R. The scale parameter of the underlying four- dimensional theory is A. In the cases that the 
Chern-Simons coupling constant takes other values, the partition functions can be also retrieved 
from the statistical model by adding another potential term for the main diagonal partitions 
in (Uni) [Si- 

The realization of partitions in terms of multiple charged partitions also plays important 
roles in multi-instanton calculus on ALE spaces [ID]; IS]- Combinatorial aspect that becomes 
key to the calculus is further elucidated in 

3.4.1 Scaling multiple charged partitions 

Let us consider charged partitions {{\^^\pr)}'rLi, where each partition A'-''-' is of order h~'^ 
and each charge Pr is of order h~^. Each charged partition is treated by the following scalings. 

a;,,(A(^)) = n(s,,|AW)/^ + C'(n°), (3.50) 

Pr = ttr/h, (3.51) 

where Sr G M>o has been introduced by v = Sr/h. 

Together with the above scalings, the correspondence leads to a scaling of the charged 
partition {fi,p) in the forms ()3.12|1 and ()3.13|1 . The density of (/i,p) is expressed in terms of 
those of {\^^\pr) as follows. 

Pbare{x\^^■,p) = j;^^Pbare{ |A^''^Pr) • (3.52) 

r=l 

We will obtain the scaling limit of this relation. We rescale a; to u by a; = Nu/h. For each 
charged partition, the inverse of u{sr\X^^^) exists, as follows from ()3.50|) . and is denoted by 
Sr{u\X^^^). It is a nondecreasing function over M, and satisfies the conditions that Sr{u\X) = u 
for u > and Sr(M|A*^^^) = for -u < — where is a certain positive constant. By using 
()3.50p and ()3.51|) . the scaling limit of each density becomes 

lim Pbare{x = u/ti lA*-^"*; Pr = ttr / K) = Pr{u — ttrlX^"^^) , (3.53) 

where 

^ rf^,^(^^ (3.54) 
du 
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As for the charged partition {fi,p), by taking account of ()3.14|) and ()3.15p . the scahng hmit of 
the density is p{u — a|/i). Therefore the relation ()3.52|) leads to 

1 ^ 

p(m - a|/i) = — ^P^(m - a^lA^*^^) , (3.55) 



r=l 



ar. 



where Na = X]r=i 

Each scaled density pr{u\X^^'^) takes values in [0,1]. We also see that dpr{u\\'^^'^) /du has a 
compact support in M and satisfies 



hoo 



,lu'-£^ = 1. (3,56) 



du 



= 0. (3.57) 

du 

At this stage it is convenient to impose the SU{N) condition ()3.46|) on the charges. This 
means that a^. are taken so that Yl^=i = 0. It follows from ()3.55|) that dp{u\p) / du is supported 
on U^^Jj., where 1^ denotes the support of dpj.{u — ar\\^''''^) / du and is thought to be a certain 
segment centered around a^. These bands become disjointed when are sufficiently separated 
from one another. In such a circumstance, by using ()3.55|) . it is possible to rewrite ()3.56p and 
()3.57p as follows. 

du^-P^ = l/AT, (3.58) 
T du 

Ir 

duu^-^ = ar/N. (3.59) 
T du 

Ir 

3.4.2 The variational problem 

We would like to consider the thermodynamic limit of each component that appears in the 
factorization ()3.47|) . In what follows, we impose the SU{N) condition on the charges. We 
also fix ttr so that they are sufficiently separated from one another. Without any loss it is 
enough to consider the case of ai < a2 < ■ ■ ■ < a^. By taking account of the previous study 
of the U{1) gauge theory, the thermodynamic limit we are investigating now is realized by 
a configuration that minimizes the energy functional with keeping a,, fixed. Therefore, the 
relevant variational problem of the energy functional should be restricted within configurations 
of the scaled densities that are expressible in terms of A^ charged partitions with the fixed 
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Figure 21: A typical behavior of the density at small h. 



charges. Such configurations may be obtained by restricting the admissible configurations 
relevantly. By observing ()3.58|) and ()3.59|1 . we require that dp{u)/du is supported on U^^Jr, 
where are certain disjointed segments, and satisfies there 

dJ4^ = l/N, (3.60) 
T du 

Ir 

duu'^^^ = aJN. (3.61) 
du 

The variational problem is formulated to find out a solution pf^*^^"* of the integral equation 
()3.27p from among the admissible configurations of the above type. 

We can also reformulate this as a problem to find out a certain analytic function: Let ^{z) 
be an analytic function that is periodic with the period l-nijR^ and behave at the infinities as 
^{z) — > ±1 as 3^2; ±oo. We suppose that $(2;) has a cut on each disjointed segments 1^ 
{r = 1, ■ ■ ■ , N) along the real axis and satisfies 

3^^$(M) = -1 on Uf_i/,, ^ ^ 

^ ^ '-^ (3.62) 

$(m) = on M \ U^^Jr ■ 

Let Cr be a contour which encircles the cut Ir anticlockwise (Figure 1^ . We further suppose 
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Figure 22: The cuts and the contours Cr on the Riemann sheet. 



that ^(z) satisfies 



R 

4-7rz 



for eacfi r. Tfien tfie solution p^^^^^ is given by the imaginary part as follows. 



dz z^(z) 



dp 



SU{N) 



dU ITX 



(3.63) 
(3.64) 



(3.65) 



3.5 SU{N) gauge theory and the Ronkin function 

We can find the solution pf^*^^"* by following the route similar to the U{1) case. Consider the 
hyperelliptic curve (j2.43|) . This curve is a double cover of C* with the branch points at two 
ends of each band Ir ()2.34p . Roots jSr of the monic polynomial Qn{x) in ()2.43p relate with 
by the conditions 

1 



2m 



zdlny = ar 



(3.66) 



Cr 



where Cr denotes the contour encircling Ir anticlockwise on the Riemann sheet. The analytic 
function 

*W = ;^(-^+ilR"'^) '^■«') 
satisfies all the above requirements including (j3.62p . ()3.63|) and (j3.64|) . Therefore we obtain 

dp^^''^\u) 1 d\ny 



du 



Nti dz ^ ^ 
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(3.68) 



The above solution can be expressed by using the Ronkin function of fsu{N)- By comparing 
(TTHHI) with (imi . we find 

pf^(^)(n) = ^ ^^^'X'^""'"^ • (3-69) 

This means 

AT^^^^^j {u, 0)=N f du iu) + const. (3.70) 

J — oo 

Let P^^^^\u) be the hmit shape of the Young diagrams obtained from by using the 

relation ()3.18j] . We can write ()3.7()j) as follows. 

Nj^^^^^{u,0) = ^(pf^W(n) +«)+ const. (3.71) 

The variational problem formulated in subsection 13.4.21 could be understood as a variant of 
that considered in |8j , where the variational problem is addressed to the dual partition function 
for the purpose of proving that the thermodynamic limit realizes the Seiberg-Witten geometry 
of the prepotential theory. However, that requires the Legendre transformation to reproduce 
the thermodynamic limit and is likely to become a detour to find a connection with the amoeba 
and the Ronkin function, while our treatment is straightforward to see the relation. 



4 Tropical geometry and crystal 

So far, we have seen that the Ronkin function of fsu{N) relates with the limit shape pf^^^^ 
as ()3.7H) under the identification ()3.(j(jj) . We expect that such a relation between the Ronkin 
functions and the gauge instanton countings is not merely a coincidence but persists further. 
To support this, we describe a certain degeneration of the amoebas and the Ronkin functions, 
and provide an interpretation from the viewpoint of statistical models. 

The parameter R is identified with the radius of the circle in the fifth dimension of the gauge 
theories. If one keeps RA fixed, one can interpret R as the inverse temperature of the statistical 
models, R = 1/T, where T denotes the temperature. This means that the large radius limit of 
the gauge theories corresponds to the low temperature limit of the statistical models. As the 
temperature approaches to zero, the statistical models get to freeze to the ground states. They 
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are determined by the charges. Crystals are complements of the ground states in the octant and 
have an interpretation as gravitational quantum foams of the corresponding local Calabi-Yau 
threefolds [33]. In this section, we argue that the low temperature limits are degenerations of 
the amoebas known as tropical geometry IM| EZI • 

The max-plus algebra or tropical semiring (R, ©, 0) is defined by 

u Q) V = max{u, v) , u&v = u + v. (4.1) 

The tropical semiring is idempotent, as follows from u®u = u. If one uses © for addition and 
© for multiplication, a tropical polynomial in two variables is defined as 

" QijU^v^ " = max(ajj + iu + jv) . (4.2) 
hi ''' 

The quotation mark in the above is used to distinguish the tropical operations from the standard 
ones. The tropical polynomials are piecewise linear functions and are responsible for some 
piecewise linear real geometry. Surprisingly, this tropical geometry can be obtained as a certain 
degeneration of the complex geometry in (C*)^. 

The idea comes from the so-called Maslov dequantization of real positive numbers. It 
is a family of semirings (M, Q)r, &r) parameterized by _R > 0. 

uQ)rV = ^Inie^"" + 6^") , uQrv = u + v. (4.3) 

For each finite R, the semiring (M, (Br, Qr) is isomorphic to the standard semiring of real 
positive number by the logarithmic map 

Log : (M>o,+,-) 0^,®r,Qr) ^^^^ 

X I — > In a; . 

On the other hand, the semiring becomes the tropical semiring at the limit i? — oo. 

u(BooV = u(Bv, uQooV = uQv. (4.5) 

By using the terminology of deformation quantization, this means that the classical semiring 
of real positive number is a quantized version of the tropical semiring. In other words, the 
Maslov dequantization says that the tropical semiring is a classical counterpart of the standard 
semiring of real positive number. 
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Anticipated from the definition of amoeba, the above dequantization deformation has a 
counterpart in amoebas. For instance, consider a Laurent polynomial f{x,y) = jbijx^y^ , 
where all the fejj's are real positive numbers and supposed to be bij = e^""^^ . As R oo, the 
corresponding amoeba Af tends to a piecewise linear curve that is described by the corner set 

4.1 Tropical limit of SU{N) amoeba 

As in the previous section, we take sufficiently separated from each other and arrange them 
in numerical order ai < ■ ■ ■ < ajy. They determine the Laurent polynomial fsu{N){x,y), given 
in ()2.28|) . by using the conditions ()3.66|) . We first study degeneration of the corresponding 
amoeba as i? ^ oo. The limit is taken with keeping ai, ■ ■ ■ ,aN and -RA fixed. This implies 
that A behaves as 0(1/ R). The thermal fiuctuations of the statistical model are suppressed at 
the low temperature, and each band of the density shrinks to the point a^. This means that, 
when R becomes very large, the parameters of the Laurent polynomial ()2.28j) are j3r ~ e^""' . 
We examine how the part of the amoeba within the strip ar < u < a^+i behaves under the 
above limit. The parts extending to the infinities can be included, by considering the cases of 
r = 0, N putting qq = —oo and aat+i = +oo. It is seen that both e^'^'x"^ (1 < s < r) and 
Q-Rasj. [r + 1 < s < N) vanish as i? — oo, if x satisfies < In |x| < a^+i. By noting this, 
we can describe effectively the zero locus which is mapped to within the strip as 

/_\N-r ^ 

y (/^A)^ ' ^ ^ 

where are chosen so that y for \y\ > 1, and y^^ for \y\ < 1. This shows that the amoeba 
within the strip < m < cir+i degenerates to two linear pieces v = ±{ru + X]fLr+i ^s) at the 
limit. 

In addition to the above, taking account of the convexity of the connected components of 
the amoeba complement, the amoeba degenerates to a piecewise linear curve that is described 
by 

N 
s=r+l 



44 



and 

N N 

[j \^{ar,v) G ; \v\ < ra, + ^ a,} . (4.8) 

r=l s=r+l 

An example of the degeneration is depicted in Figure l23l 




Figure 23: Tropical degeneration of the the amoeba of fsu{2) ■ 

We turn to consider the behavior of the Ronkin function at the tropical limit. Over the 
amoeba complement, the Ronkin function is the piecewise linear function Sf^^^^y At the 
tropical limit, where the amoeba degenerates to the above curve, these two functions coincide 
all over M^. Let us write the piecewise linear function thus obtained as 

Soo{u,v) = lim Nf {u,v). (4.9) 

We can reproduce this 5*00 from the piecewise linear curve ()4.7|) and ()4.8|) . except an addition 
of a constant to the function. This is because the piecewise linear curve is the corner set of 
5*00 and the gradient VS'oo equals to the corresponding lattice point of the Newton polygon on 
each connected component divided by the curve. The constant may be found, for instance, by 
considering the Ronkin function over E^^ where a is a vertex of the Newton polygon. In this 
manner we obtain 

Soo{u.,v)= max ( -rf^ + ((«, "w), «) ) , (4-10) 
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where 

d(ifl)=Y!r=l(^r {l<i<N), ^^^^^ 

d{o,o) = c?{o,±i) = . 

The set of {u, v, Soo{u, v)) becomes a facet of a three-dimensional polyhedron. We associate 
a three-dimensional lattice element with each integer lattice point of the Newton polygon by 

v^ = i-aA) eZ\ (4.12) 

where a G Af^^^^^^ fl Z^. By using these Va, we can describe the polyhedron as follows. 

Psu(N) = {iu,v,w) eW"; {iu,v,w),v^) + d^>0 Va G A^^^^^, n } , (4.13) 

where ( , ) denotes the standard inner product on R'^. By comparing the above with (j4.1LH) . we 
see that the facet of this Vsu(n) is actually given by 5*00. See Figure EH 




Figure 24: The three-dimensional polyhedron Vsu{n)- 
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4.2 Tropical geometry and quantum foam of local geometry 

The polyhedron Vsu{n) has an interpretation in the local geometry and accords with the quan- 
tum foam picture [IH],|niI- We introduce three-dimensional cones do, cii, ■ ■ ■ , a2N-i by 

CT2i = ]R>ot^(i,0) + K>0t'(i+1,0) + IR>0^^{0 -1) , 
0"2i+l = IR>0'"(i,0) + IR>0'^(i+l,0) + IR>0'W(0,1) • 

All these cones and their faces constitute the fan that describes the local SU{N) geometry X. 
Two-dimensional cones of the fan determine two cycles of the geometry as closed subvarieties 
that are invariant under the torus action. Cones that are generated by f (fc,o) and "^(0,-1) deter- 
mine vanishing cycles in the fibred ALE space, where 1 < k < N — 1. Cone generated by f (1,0) 
and f(o,o) determines the base P^. 

The polyhedron ()4.13p . more precisely, that rescaled by 1/h, emerges naturally when quan- 
tizations of the local geometry X equipped with a Kahler two form uj on it, are considered. 
For the consistent quantizations, this u needs to be quantized. Topological A-model strings 
suggest the rule that u is quantized in the unit of string coupling constant gst |TH] . 

— [u]eH\X,Z). (4.15) 
Qst 

This means that the Kahler volumes of two cycles become integral. 

Tb = — [ CO eZ>o, Tk = — [ to eZ>o (A: = 1, ■ ■ ■ , AT - 1) , (4.16) 
9 St Jvb 9 St Jvk 

where Vf, and Vk denote respectively the base and the vanishing cycles in the fibre that 
correspond to the cones generated by f(fc,o) and '^(o,-!)- Relation between quantizations of the 
local geometry and the statistical models is found in [3^. In particular, the integral parameters 
()4.1(ij] are converted to 

9st = Rh, n = -—log{RAY^ , Tk=Pk+i-Pk, (4.17) 

9st 

where Pr are integers arranged in numerical order Pi < P2 ^ ■ ■ ■ Pn- We will associate these 
integers with partitions subsequently. We will also impose the SU{N) condition on the integers 
and slightly restrict allowed values of the parameters so that Ylk=i ^-^fe ^ A^Z>o. 

When the limit i? ^ 00 is taken with keeping RA fixed, it makes Th vanish under the 
identification ()4.17j) . This implies that we can put = from the beginning, for the present 
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purpose. Let uj be such a quantized Kahler form. The geometric quantization of (X, uj) requires, 
first of all, a holomorphic line bundle L on X which the first Chern class equals ci(L) = [uj]/ Qst- 
We may take such a line bundle as L = 0{D), where D is a certain toric divisor. Rays or edges 
of the fan determine closed subvarieties of codimension one that are invariant under the torus 
action. The divisor is a formal sum of these subvarieties with integral coefficients. 



D = J2PaDo., (4.18) 



where is the subvariety that corresponds to the ray generated by v^- The integral parameters 
()4.16p relate with the above Pa by using ci{0{D)) = [ouj/gst. The identification ()4.17p allows 
us to express the coefficients pa by the charges as follows. 

P{ifi) = El=lPr (l<^<Ar), ^^^^^ 

P(o,o) = P(o,±i) = 

Physical states in the geometric quantization appear as the global sections of 0{D). Thanks 
to that X is a toric variety, these states are labeled by integer lattice points in a convex 
polyhedron determined by D 



V{X,D) = ^^{u,v,w) eM.^ ] {{u,v,w),Vc,) +Pa>0 Va}. (4.20) 

Thus space of the physical states has a basis Xm, where m runs over V{X,D) fl Z^. Each 
Xm is interpreted as a quantum of {X,u). The charges Pr are scaled to by pr = ar/h at 
the thermodynamic limit or the semiclassical limit. Comparing ()4.19p with ()4.1H) . we find 
Pa = da/h. This means that two polyhedra V{X, D) and Vsu(n) are similar. They relate with 
each other by the similarity transformation. 

V{X,D) = ^VsuiN). (4.21) 

The number of the physical states is infinite since the cardinality of V{X, D) flZ^ is infinite. 
In order to manipulate such an infinity, another polyhedron that originates from a local singular 
geometry is introduced [SI]. 

Vs^ng = {{u,v,w)eW'; {{u,v,w),v^)>0 « = (0, ±1) , (X, 0) , (0, 0) . } . (4.22) 
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The underlying singular geometry is a Za? orbifold of (9 © 0{—2) 
fan consisting of three-dimensional cones ri,r2 and their faces. 

Ti = R>o'y(o,0) + I^>0'y(Af,0) + I^>0'W(0 -1) , 
T2 = M>ot;(o,0) + IR>0'i^(Ar,0) + l^>0'i^(0,l) • 




Figure 25: The three-dimensional polyhedron Vsing- 

The above polyhedron includes V{X,D) as an unbounded subset. In [Hlj, the number of 
the physical states is regularized and is prescribed to be the cardinality of V^{X, D) \~\ I?, 
where V^i^X^D) denotes the complement Vsing \ 'P{X,D). Since it becomes bounded, this 
relative counting gives a finite answer. When the charges become very large, the cardinality is 
approximated by the volume measured by the volume element dudvdw. Therefore, the number 



P and is described by the 



(4.23) 
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of the physical states becomes ~ Vol(P'^(X, D)). The volume can be computed as 

1 ^ N ^ 

Vol{V^{X, D)) = - Y^iPr - Ps)' + jJ2Pr^ (4-24) 

r>s r=l 

which turns to equal —D^/3\, where denotes the self- intersection number obtained by in- 
tersecting D with itself three times. 

The counterpart of the singular geometry is a Laurent polynomial consisting of monomials 
1, x^, y, and reflecting that the geometry is the Z^v orbifold of C © 0{—2) ^ As such, 
we will take 

fsTiN)i^,y) = - 1 - (i?A)'^ - iRAf{y + y~') . (4.25) 

The amoeba of fsuiN) tentacles which asymptote to the straight lines ()2.3H) but there 

appears no hole. The tropical limit of the Ronkin function Nfbase becomes 

■IsUiN) 

5^-(u,t;)= max ({{u,v),a)). (4.26) 

a=(0,0),(Af.O),(0,±l) V / 

This piecewise linear function describes the facet of Vsmg- The counterpart of V^{X, D) is the 
complement '^^[/(jv) = 'T^sing \ 'Psu{n)- This is similar to V^i^X, D). 

V%X,D) = ^V^suiN)- (4.27) 



4.3 Tropical geometry and crystal 

By using the correspondence ()3.44|) . charged empty partitions describe a specific partition, 
which is called core. Let us consider A^ charged empty partitions {(0,Pr)}^i that charges 
are of order and satisfy the SU{N) condition ()3.4fjj) . The corresponding A^ core is denoted 
by 

f^N-core- See Figure |27| Asymptotics of f^N-coie &t the limit h — > is obtained by scaling 
each charged empty partition according to ()3.50p and ()3.5H) . In particular, the charges are 
scaled by Pr = ar/h. The scaled density of fiN-coie becomes 

1 ^ 

PN-coTe{u) = p{u\nN-core) = ~ ^r) ■ (4.28) 

r=l 

50 




Figure 26: The complement P|[7(jv) = 'Psing \ VsuiN)- 

Each charged empty partition contributes to the above as a step function. Comparing ()4.28|) 
with ()4.10|) . one finds that this pN-cove is written by using the piecewise hnear function 5*00 as 

1 dS^ju, 0) 

Piv-core(^) = ^ ■ (4.29) 



The above expression suggests a certain relation between the core and the tropical limit 
of the SU{N) amoeba. To see this, we first compute the energy of fiN-coie- The charges are 
arranged pi < ■ ■ ■ < Pn- We evaluate the energy functional (j3.23|) at p = pN-coie- The function 
■y{u;A;R), that is defined by the conditions ()3.24|) and ()3.25p . has the expression 

7(«; A; R) = ^u' - log(/?A)«2 -^""^^ ^^^^^ ~ Lz,{e'^^)) . (4.30) 
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Pn - Pn-1 



P2-V1 



Figure 27: /xtv- 



By using this, the energy of fiN-core becomes 



E[p 



R 



N 



N 



N 



Af-corcJ 



- <^ ^(a, - asf + AT ^ an _ iog(i?A) ^(a, - 

(_ r>s r=l J r>s 

J2(^r - as) + E(C(3) - ^^3(e-^(-"-))) . (4.31) 



6^ 



r>s r>s 

This is the perturbative prepotential of five dimensional = 1 supersymmetric SU{N) Yang- 
Mills on X 5^. The first term of (j4.3ip dominates when the limit i? — 00 is taken. Therefore, 
taking account of ()4.24|) and ()4.27|) . we obtain the following estimate at the low temperature 
limit. 



E[p 



N-core\ 



R-Yo\(v% 



(4.32) 



A set A^(pi,- ■■ denotes a set of plane partitions that main diagonal partitions are 
fJ^N-core- It follows from ()3.19j) that the energy ()4.31|1 is the free energy of random plane 
partitions restricted within A^(pi, ■ ■ ■ ,Pn)- 



(4.33) 
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The ground state of the above restricted model is a plane partition that minimizes |7r| in 
A^(pi, ■ ■ ■ ,pn)- Such a plane partition is determined uniquely by the charges and is denoted 
by ttgpp- Elements of this ncpp are expressed by using the core as follows. 

'KGPp{m)i = UN-core i+\m\ ■ (4.34) 

The number of the cubes of tcgpp becomes 

Et 2 
/^TV-core i 

i=l 

= jT.(P^-Psf + ^T.Pr + Oip')- (4.35) 

r>s r=l 

At the limit ^ — > 0, by taking account of ()4.24p and ()4.27|) . we obtain 

lim IttgppI = N ■ Vo\(v'su(N)) ■ (4.36) 

The plane partition t^gpp becomes the ground state of the original random plane partitions. 
At the low temperature limit the statistical model freezes to hgpp- There is a description j^l] 
that the complement Vg^^j^-^ is identified with the plane partition ttgpp at the thermodynamic 
limit. Consider the linear transformation A in M^. 

/ -1 l\ 



A 



(4.37) 



1 1 
\ -N 1 / 

The image of the polyhedron Vsmg is the positive octant O = {{u,v,w) G M^; u,v,w > 0}. 
Thus the image of the complement Vg^^j^^ is 

A{V'su(N)) = 0\AiVsuiN)). (4.38) 

This AiVgm^j^^) can be identified with the shape of the plane partition tcgpp- We assemble a 
plane partition vr in the positive octant O, by a rule slightly different from the standard one. We 
substitute unit cubes by rectangular solids of the size \/2 x \/2 x 1, and number 7r/4-rotated 
squares in the first quadrant of the {u, t')-plane as in Figure |2H1 We stack pieces of the 
rectangular solid vertically on each (z, j)-square of the first quadrant. As the result, the plane 
partition itgpp is assembled in the positive octant O. It is straightforward to see that this itgpp 
scales to A{Vgjji^j^^) at the limit h ^ 0. 

53 




Figure 28: tt / i-rotated squares in the first quadrant of the {u,v)-plane. 
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